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ABSTRACT 


In this study, the governing differential equations 
and the corresponding finite difference equations, for one 
dimensional and two dimensional horizontal flow models, 
have been developed to determine numerically the steady 
state velocit37' field and the resulting temperature distri- 
bution in the water body, which receives the heat rejected 
from a thermal source and dissipates this energy to the 
environment. Tiie flow pattern for the two limiting cases 
of inviscid and creeping f3.ow, as well as for flow at 
finite low Reynolds numbers, have been analysed. The 
effect of various parameters, such as the aspect ratio of 
the water body, intake and outfall locations, flow rate, 
eddy diffusivity etc., have been investigated. The study 
provides guidelines for limiting the power plant capacity 
for a given water body and determines the effect of the 
power plant capacity on the intake temperature rise. 



CliAPTER I 


INTSODU'CTION 


1.1 The Problem 

The importance of heat rejection to water bodies 
is primarily due to its use in thermal and nuclear power 
plants, where large quantities of cooling ■vj'ater are utilised 
to condense the esdiaust steam from the turbine as water. 

Tor an efficient and satisfactory operation of these plants, 
it is essential to have an efficient system for heat 
re j ection. 

The continued growth of electric power supply is 
vital for our industry and, thus, for our national economy. 

It is imperative that this growth be accomplished economi- 
cally and v/ith an environmental impact idiich is on accept- 
able level. As the growth of industry based on electric 
power continues, unit sizes, plant capacity and site locations 
present new challenges. The bifLk of the electric power, at 
present, is being generated through the use of steam driven 
thermal povrer plants and these p3.ants use large quantities 
of cooling water to condense the exhaust steam from the 
turbine. This keeps the temperature of condensation low 
so as to maximize the energy conversion. The condensate 
is recirculated back to the boiler and the cooling water 
is generally heated by 8 to 12 C, through energy exchange 
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in the condensers. One of the most important factors to 
he considered in selecting plant locations is the method 
of waste heat rejection and the availihility of water 
resources to accomplish it. 

Cooling water is generally obtained from natural, 
or artificial, lakes or from rivers, and is returned to 
the source at a higher temperature. Ultimately, almost 
all this energy from the cooling water, is dissipated to 
the atmosphere through energy exchange at the surface, with 
a negligible loss to the bed of the lake or the river. 

This rejection of heat from the cooling water alters the 
naturally existing temperature balancd in a body of water. 

The heat transfer mechanisms w'dch occur at the 
surface of a cooling water bod37- comprise of the following? 

(a) Incoming short-wave solar radiation, 

(b) incoming long-wave atmospheric radiation, 

(c) portions of both short-wave and long -wave radia- 
tion which are reflected by the water surface, 

(d) long -wave back radiation from the water to the 
atmosphere, 

(e) heat exchange due to conduction and convection, 

(f) heat loss due to evaporation and gain due to 
condensation . 
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The comhined effect of all these heat transfer 
mechanisms determines the rats of cooling for water hodies 
under specific conditions. The last three of the above 
heat transfer mechanisms, i.e., back radiation, conduction- 
convection and evaporation are functions of the water body 
surface temperature. This water temperature is subject.-d 
to control through variations in plant design or operating 
practices. Although the quantity of heat to be rejected 
is relatively fixed for a certain plant c8.pacity, condenser 
design and water flow rate can be manipulated to achieve 
different values of discharge temperature and flow velocity. 

The maximum temperature increase at any given point 
and the affected area or volume can bo estimated for speci- 
fic conditions. The maximum ternpernture can bo depressed 
by mixing the returning cooling water as thoroughly as 
possible with ambient water. At the other extreme, the 
cooling v/ater may be floated over the ambient water as a 
stratified layer, maintaining high surfa.ce temperature. 

The latter scheme has the advantage over the former in 
that the rate of heat loss to the atmosphere is maximum, 
and, therefore, tLie area affected is minimum, even though 
local temperatures are higher. At the same time, the 
intake for the power plant has to be located on a lake or 
idver, so as to obtain water at the lowest possible tempe- 
rature for a high thermal efficiencjr and a discharge 
located to pre'B'ent warm water recirculation to the intake. 



1.2 Previous Work 


Not much work has so far "been done in this area. 
Initially, this area was considered mainly due to the 
thermal pollution effects arising from heat rejection. 

The need of a fluid flow and heat transfer analysis has 
been realised only during the last few years, with a view 
to improve the power plant efficiency. 

Steady laminar flow, with closed streamlines at 
large Re 3 molds number, was analyzed by Batchelor [!I‘ The 
effects of various factors on the amount of recirculation 
of water, after it has been used as a coolant in the thermo- 
electric plants and returned to the river, or lake, from 
where it was originally drawn, vras perhaps first carried 
out by Bata L2D . yome of the more important elements 
which need be considered in a well designed circulating 
water system were highlighted by Richards D3. 

Raphael Zhl presented a procedure for predicting 
the temperature of various water bodies, such as shallow 
lakes, flowing streams and detention reservoirs, from 
weather records, the inflow and outflo\; characteristics, 
the surface area and the volume of the water body. This 
procedure assumes that the thermocline is absent and the 
water so stirred by wind or internal currents that the 
temperature is uniform throughout. A mathematical model 
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for predicting temperatures in rivers and reservoirs 
was developed by Delay and Seaders C53 They computed 
the net monthly energy exchange from the reservoir and 
distributed it vertically so that the resilLting tempe- 
rature gradients approximated the standard for each 
month. Dake and Harleman Z7l formulated a theoretical 
model for the time dependent vertical temperature dis- 
tribution in a deep lake during the yearly cycle by 
using heat-flux balance at the water surface, which 
accounts for back radiation and evaporative loss, as a 
boundary condition. 

It was realised later on that to determine the 
methods for minimizing the recirculation of vraste heat, 
in the design of cooling water systems for thermal power 
plants, it was desirable to understand the nature of 
mixing of the -^^aim water discharge. Accordingly several 
experiments on the mixing of the buoyant jet of water 
being discharged horizontally at the surface of a body of 
water were performed in tho laboratoiy by Yuan, VJiegel 
and Ismail and by Tamai, Wiegel and 
The results obtained from these experiments compared 
favourably with the actual data obtained from different 
thermal power plants. Subsequently it was shown by 
Stefan 0^* '1311 "tbat the flow and cooling of a warn 
water surface layer could be represented in a physical 
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model if, first, flow of water in and out of the model, 
and second, atmospheric conditions, it was exposed to, 
were controllable. He simiiLated the phenomena of mixing 
at the outlet, stratification, and the surface heat 
transfer through the use of physical models, in experi- 
mental design to simulate a heated water discharge from a 
channel into a deep lake for lateral spread was also 
developed. 

Brown |33 discussed the various methods for the 
disposal of waste heat and concluded that studies were 
needed to develop improved and less costly cooling device 
that would be required increasingly at new plant locations. 
He also emphasised the need to explore the possibilities 
of a beneficial use of the waste heat. Uie alternative 
methods of heat dissipation from power plants and the 
economic evaluation of these alternatives was presented by 
Hauser Loius Analysis for reducing the cooling pond 

area and wa.ter loss, by designing the pond for higher opera- 
ting temperatures, was presented by Brue and Alden m- 

The thermal effects of povrer plants on lakes in 
terms of its temperature cycle were studied by Moore and 
Jaluria 153 ^ one dimensional model based on the 

assumption of a constant Bicha^rdson number at the base of 
the stratified layer was developed. The model gave a 
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constant heat flux across the the mio dine. The model was 
then perturbed in terms of heat flux as well as the 
thermal diffusivity bo give the power plant impact on a 
lake. Based upon the assumption of horizontal isotheims 
at all times, Huber, Earleman and R 3 ran |j6 3<3.eveloped a 
mathematical model to predict the vertical temperature 
distribution in stratified reseivoirs. The model included 
the effects of distribution of heat dtiiin the reservoir, 
by convection and diffusion heat sources and sinks at 
boundaries, and internal dissipation of solar radiation. 
Temperature profiles predicted from the solution of this 
mathematical model agreed well -ith measured values in both 
the field and the laboratory, Hindi ey and Miner 
observed that the temperature distribution in a cooling 
water body, being employed for a power plant, was dependent 
upon the manner in w.hich the heated effluent mixed with the 
receiving water body and upon the rate of heat exchange 
between the water surface and the atmosphere* It was found 
that, far from the receiving end, where mixing was nearly 
complete, the dominant cooling mechanism was heat exchange 
with the atmosphere. The rate of this cooling was directly 
proportional to t]ie surface heat exchange coefficient, 
which in turn is a function of several environmental varia- 
bles, particularly the wind velocity. 
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Grubert and Abbott QlS]] foimulated a mathematical 
model for nearly horizontal stratified flows. They showed 
that the characteristic directions in a stratified fluid 
divide into pairs, each pair being associated with a fluid 
layer. Miyazaki D9J studied the heated two dimensional 
jet discharged at the water surface and carried out an 
analysis at various Eichardson numbers for determining the 
velocity and temperature distributions. 

Jaluria, Variyar and Mehta {^20, 21]3 formulated a 
mathematical model to predict the temperature and velocity 
profiles in a body of water, due to the basic mechanisms 
in the presence of wind, cloudiness, back radiation etc. 

They computed the equilibrium average surface temperature, 
which is defined as the temperature attained by the surface, 
if the ambient conditions are held constant at specific 
values; and studied the heat transfer mechanisms at the 
surface in detail. They also developed a more generalised 
model, considering convection effects caused by the outflow 
and inflow of cooling water and the heat addition from the 
power plant. Certain specific, simplified, one and two 
dimensional models were studied by them in detail. 

1,3 Present Work 

In the present work, one dimensional and two 
dimensional, horizontal flow models have been considered. 
Corresponding finite difference equations have been 
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developed to determine numerically the steady state velo'city 
field and the resulting temperature distribution in the 
^'jater body, which receives the heat rejected from a thermal 
source and dissipates this energy to the environment. These 
models are related mainly to the problem of heat rejection 
from a power plant to a water bodiy, such as a lake, river, 
sea etc. These may also represent the problems of heat 
rejection from the condensers of an airconditioning plant to 
a water pond or the heat rejection from industrial systems 
to water bodies. 

A nonlinear fourth order partial differential 
equation, in the form of nondlmensionalized variables, is 
obtained for the flow, in a tiro dimensional model for the 
water body, neglecting varia,tions in the vertical direction. 
The finite difference formulation of this equation in the 
explicit form was used for computing the recirculation at 
different Rejmolds numbers. The recirculation was also 
computed for the two limiting cases, of inviscid and creeping 
flow. The effect of various parameters, such as aspect 
ratio of the water body, the intake and outfall locations 
and eddy viscosity, on the nature of recirculation was 
studied. The resulting velocity field was obtained and 
the general nature of flow discussed in terms of stream- 
lines and velocity distributions. 
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For finding the temperature distrihution, a one- 
dimensional model, neglecting variations in the vertical 
and lateral directions, was initially studied. The effect 
of heat rejection from the power plant on the temperature 
distribution in the water body was computed and its depen- 
dence on various parameters, such as, length of the water 
body, temperature difference between intake and outfall 
temperatures, flow rate etc,, have been investigated. The 
effect of turbulence has been considered in terms of the 
eddy dlffusivity The net surface heat exchange was 

computed by using a method similar to that outlined by 
Jaluria, Variyar and Mehta [1203, vrith respect to ambient 
conditions such as solar flux, wind speed, ambient tempe- 
rature, cloudiness factor, etc. 

In the second phase of the present work, a two- 
dimensional model, neglecting variations in the vertical 
direction, was developed to determine the temperature 
distribution in the water body. The temperature distribu- 
tion has been computed for creeping flow and for flows at 
small values of the Reynolds number. The impact of the 
intake and outfall locations on recirculation was studied 
so as to obtain a system with minimum intake temperature 
rise, for rejecting a particular quantity of theimal 
energy. The effect of various parameters, such as dimen- 
sions of the water body, intake and outfall temperature 
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dd-fferencfi^ floar rate, etc,, on the temperature distribution 
has been considered in detail* The physical significance 
of these results with respect to actual systems is discussed. 

All the computations were carried out on IM 7041f- 
at I,I,T, Kanpur. The meteorological data needed for the 
study has been obtained largely from Indian Meteorological 
Department and partly from the analysis given by 
Raphael 



CHAPTER 2 


FOIMTJLATION OF THEOEETICAL MODELS AND METHOD OF SOLUTION 

The behaviour of a water body, which acts as a 
cooling water source for a power plant, is determined, to 
a large extent, by the nature of flow field set up due to 
the outfall of warn water from the condensers of the power 
plant and by its interaction with the environment, to 
which it ultimately dissipates the heat rejected by the 
power plant. Hence, in developing an analytical model for 
a cooling water body, the mathematical representation of 
the various mechanisms underlying its energy exchange with 
the environment is of considerable importance. 

This chapter considers the analysis and the methods 
of solution for different cases of fluid flow in the water 
body, under steady state conditions. The basic features and 
assumptions made in carr 3 ring out the analysis are outlined 
as follows? 

(a) Finite dimensions of the vater body. 

(b) Energy transfer takes place at the intake and 
outfall of the cooling water for the power plant 
and at the surface of the water body. Energy 
transfer from the sides and the bottom of the 
water body is neglected. This assumption is 
adequately supported by various investigations. 



13 


(c) An idealized configuration, of vertical sides 
and flat horizontal "base and top surface, is 
considered. This can eventually be modified to 
consider, as close as possible, the exact topo- 
logy of a given water body. 

(d) Though ambient conditions are functions of time, 
they are assumed to have specific values over a 
considerable period of time. 

(e) The effect of viscous dissipation is neglected. 

(f) The physical properties of the fluid are assumed 
constant, with respect to temperature variation. 

Any other approximations employed in addition to the above, 
have been mentioned in the course of analysis, 

2,1 Recirculation Models 

The amount of recirculation, to the cold water 
intake, of the heated water directly affects the efficiency 
of operation of a power plant and also usually lowers its 
maximum capacity. It depends on the geometry of the 
intake, the design of the outfall, the distance between 
the intake and outfall, the discharge and mean depth of 
the flow in the water body, the relative amount of diver- 
sion and the degree of heating of the diverted water. 
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The analysis proceeds "by applying the principle 
of conservation of mass, momentum and energy to a diffe- 
rential control volume at any point in the water body. The 
basic equations, for incompressible floi-/ under steady state 
conditions, ares 

= 0 ( 2 . 1 ) 

( V. V ) V = 1 ( F - Vp) + V V (2.2) 

( V.V ) t = Vjj t (2.3) 

where V denotes the velocity, F the body force, t the 
tempera turej v the eddy viscosity and the eddy diffu- 
sivity. The net amount of heat exchange per unit area with 
the environment appears as a boundary condition at 

the surface. 

The solution of the above equations gives the 
velocity field and the temperature distribution in the 
water body. The equations are quite complex and are usually 
coupled through the temperature t because of the buoyancy 
force which arises from the ( F - Vp) term and is a func- 
tion of t. Thus, these eqtiations need to be further 
simplified before a solution is attempted. 

A two dimensional model with complete mixing in 
the vertical direction is considered. This situation 
frequently arises in the case of a shallow water body. 



The intake and outfall are assumed to be located at the 
same level, since due to vertical mixing a vertical separa- 
tion can not be considered. The vertical direction being 
fully mixed, the buoyancy effect of the temperature varia- 
tion with depth does not arise. As a result, the problem 
becomes a decoupled one, i.e., the flow pattern may be 
determined without considering the tem^Derature field. The 
temperature variation in the horizontal plane may then be 
determined from the known flow field. The basic equations 
become s 
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where Q is the net amount of heat exchange per unit 
surface area and 11 the depth of the water body. 


2,1.1 Governing Equations 

The eddy viscosity at any point in the flow 
region for turbulent mixing in the directions x and y 
is assumed to be equal so that 


V 



16 


The equation (2.4) can he satisfied identically 
by using a stream function given by 
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9 il/ 
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( 2 . 8 ) 
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Tho values of u and v in terms of the stream 
function are substituted in Equations (2.5) and (2,6). Then 
to eliminate the pressure term pj equation (2.5) is diffe- 
rentiated v±th respect to y and equation^ (2.6) with 
respect to x and the resulting equa^tions are subtracted 
from each other to yield the follo\- 7 ing equation - . 
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The equation (2.10) is non-dim ensionali zed by 
using the following dimensionless parameters s 
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where U is the uniform velocity at inlet to the water 



17 


body and 1 is the opening at the inlet point. Here P 
represents the dimensionless stream function. However for 
convenience of notation, this has been shoim on Figures 
through ^ , 

The non- dimen sionalized equation thus obtained is 
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3 Y ' a X 
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(2.13) 

where Re = ESyxiolds nunber = ~~ (2.14-) 

The solution of Equation (2.13) alongwith the 
requisite boundai^;^ conditions gives the velocity field for 
viscous flow. The equation is a nonlinear, partial diffe- 
rential equation and its general solution is not known. It 
is amenable to a numerical solution in finite difference 
form, being somewhat simpler at low Re3molds number. The 
solution can be obtained for three different cases outlined 
below. 


(a) Inviscid Flow 

The stream function for two dimensional inviscid 
flow, being irrotational, satisfies the Laplace equation 

IL • 0 • 

’ = 0 (2.15) 

with the following boundary conditions, as shown in 
Fig. 1. 
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This implies constant value at the walls as shorn in 
Fig. 1 . 
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(b) Creeping Flow: Re << 1 

If the Reynolds number is -veiry small, the R.H.S. 
of Equation (2.13) dominates the flow and the inertia 
terms on the L.H.S, can be neglected, lliis is the special 
case of creeping flow, where viscous effects predominate 
and inertia terms are negligible. The equation then 
reduces to 

4 _ 

V ?= 0 (2.16) 


This is a fourth order linear equation Imown as 
the biharmonic equation in two dimensions. The solution 
of this equation alongwith the follomng boundary condi- 
tions gives the velocity field for creeping flow. 
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and = Lj 0 < x < W ; = -r^ = 0 

— o A 0 X 

This implies constant value of Tp at the walls and also the 
value of the gradient as given above. 

(c) Flow at Higher Reynolds Numbers 

For obt3-ining the velocitj^ field at differentx 
Reynolds numbers. Equation (2.13)} with all the terms, is 
solved, alongwith the boundai^r conditions given by equation 
(2.16 a). Even for very large Reynolds numbers, the R.H.S 
of Equation (2.13) can never be negligible near a solid 

k _ 

boundary because the no slip condition forces V ip to be 
very large, of order Re, near the wall, 

2.1,2 Method of Solution 

Equations (2.13), (2.l5) and (2.16) are solved by 
using finite difference techniques. Employing a rectangular 
mesh, the surface area of the water body is divided into M 
parts along the y ~ axis and N parts along the x - axis, 
each part being of length ax. The various partial 
derivatives are replaced with their central difference 
approximations, which convert the governing differential 
equations into a set of simultaneous algebraic equations. 
Appendix I gives the finite difference form of the equa- 
tions. The finite difference equations are solved by 
using the Gauss - Siedel iterative method. 
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Before starting the iteration sehemej all the mesh 
points need to be initialized. The computer time needed for 
achieving a given criterion for convergence depends on this 
initialization and as such the guess for this should be as 
close to the exact solution as leossible. 

The rate of convergence for inviscid flow is much 
faster than for other flows due to the lower order of the 
governing equation. Therefore, the Equation (2.15) -3 
is solved first by using the following initialization. 

0 

= 0 for all (i, j) 

^ j <] 

where (i, 3) represents the interior mesh points in the 
flow area. 

The solution obtained for inviscid flow is used 
as an initial guess for solving the equation (2.16). 

Finally the creeping flow solution is used as an initial 
guess for solving Equation (2.13)* The solution of 
Equation (2,13) Tor higher Reynold numbers is obtained by 
using the value of stream function for lower values of 
Reynolds numbers, close to the new values, as an initial 
guess. 

__0 

Starting from these initial values of 1 ’ 

1 2 r ’ 

the successive iterates ^i,j» ^i,i’ ’ 

are obtained by using the finite difference equations 
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given in Appendix I. The iteration is continued till 

convergence is achieved i.e., the difference between 
_r+1 _r 

- . and 4' . . becomes smaller than a specified 
value for all values of i, o, where r is a given itera- 
tion. 


2.2 Energy Exchange at the Surface 

The energy exchange at the surface is due to 
various heat transfer mechanisms outlined below s 

(a) Eadlation from the Sun and from the -Atmosphere 

At a point outside the atmosphere of the earth, 
radiant energy is received from the sun, on a surface 
normal to its rays, at a rate of approximately 1165 KsCpl/ 
hr/Sq. m Tills value is called the 'Solar Constant', 

A plane normal to the rays of the sun at the surface of 
the earth receives considerably less energy as compared to 
the solar constant because moisture, gases and solid par- 
ticles in the air scatter and absorb much of the incident 
radiation. Although much of tlie radiant energ37' is lost, 
some of this energy is recovered as diffuse radiation. 

Of the total quantity of the direct, or short 
wave, radiation that reachefj tlie surface of a water body, 
a portion is returned unchanged due to reflection at the 
surface and scattering by bubbles and suspended particles 



immediately below it However, for practical engineer- 

ing computations, solar radiation and reflected radiation 
can be combined in one function, (Qj_ - or s 

Effective absorbed solar radiation = Q.. = Q^- - Q„ 

o J. J. 

If c is the average cloud cover in tenths of 
sky covered, tho net short wave radiation is obtained from 
the correlation given by Raphael 

Qg = (1 - 0.0071 c^) - Q^) (2.17) 

(b) Back Radiation 

Effective back radiation m.3.y be defined as the 
difference between Q , the longwave radj.s.tion leaving a 
body of water, and Qg^, the longwave radiation from the 
atmosphere being absorbed by the body of water Q+Q. Atmos- 
pheric radiation does not follow any simple law as it is a 
function of many variables, such as the distribution of 
moisture, temperature, ozone, carbon dioxide, etc. Thus, 
this has been found to be dependent on cloud height, as 
high," middle and low clouds, and the amount of cloud, as 
scattered cloud, broken cloud and overcast. However, as 
most weather observations give cloud amount in tenths of 
sky obscured, the variation of the atmospheric radiation 
factor has been given as a function of cloud amount and 
vapour pressure by Raphael [30 • effective back 
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radiation, can be computed as, 

V = - «a = 0-97° ^ (S-IS) 

•where a is Stefan-Boltzman constant, is absolute tempe- 
rature of water surface, Tg is absolute ambient temperature 
and e is sky radiation factor, tabulated by Paipliael ra. 

(c) Convection 

Sensible heat is convected to or from a body of 
water whenever a temperature difference exists between air 
and water. The basic equation for the convection heat 
transfer, Q^, from the water surface to the air is, 

where h is the surface heat transfer coefficient. Hie 
value of h is largely dependent on the condition of the 
air in contact x/ith 'the water and hence it is a function 
of the wind velocity, its temperature, direction and 
turbulence level etc. The following correlation by Wada 
[2^ has been found to be most s’Uttable for conditions 
relevant to India, 

h = 2.77 x 10~^ (0.48 + 0.272 V ) K Cal/m^ sec °C 

( 2 . 20 ) 

where V,, is the wind velocity along the water surface 
in me ter/ sec. 
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(d.) Evaporation 

Vtoen tile vapour pressure of ambient air is less 
than the saturated vapour pressure, at the vmter surface 
temperature, water evaporates into the air, removing heat 
from the water mainly due bo the energy required to evapo- 
rate the water and to a sm.all extent by the secondary heat 
transfer due to the sensible heat contained in the water 
removed by evaisoration. Tiie heat transfer by evaporation 
is given by Hindley and Miner D7D 

Qg = 2h [> ® ^^a^J (2.21) 

where e (Tg) is the partial vapour pressure at water sur- 
face temperature in mm of Hg, e (Tg^) is the partial vapour 
pressure at ambient temperature in mm of Hg and h is 
the heat transfer coefficient given b 3 /' equation (2.20), 
Also 

e (Tg_) = Gg (T^) X R.E. 

where eg (Tg_) is the saturated vapour pressure at ambient 
temperature and R.H, is the relative humidity at ambient 
temperature. 

The rate of net heat loss, Q, from the surface 
is given by 

Q = Qg + "^c (2.22) 
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I'Jhen heat is being rejected by a power plant to a 
water body, the average surface temperature rises and the 
equilibrium condition is obtained i/hen 

% ~ *^e %r “ ^s ^c ^ ~ % (2,23) 

where Qp is the rate 01 heat rejected 'bj the power plant 
per unit area of the water body surface. 

2.3 One Dimensional Model 

Physically, a one d±mensional model represents 
the circumstance inhere mixing in the water body in two 
directions is so vigorous that excha.nge of heat and mass 
need be considered only in the third direction. One 
dimensional horizontal model with energy exchange at the 
surface and variation in the longitudinal direction, is 
considered. Therefore, temperature in the vertical and 
lateral directions is assumed to be uniform. The theore- 
tical model, under steady state conditions, is developed 
below. 


2.3.1 Gove rning Equa ti on s 

The surface of the water body is assumed to 
consist of a number of finite divisions, each of width 
&x. The finite difference equations are obtained by 
considering the energy balance for each element i.e., 


Ehergy input 


Energy output 


Losses 
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(d) Evaporation 

VJhen t]iQ vapour pressure of amtient air is less 
than the saturated vapour pressure, at the x^rater surface 
temperature, -wator evaporates into the air, removing heat 
from the water mainly due bo the energy required to evapo- 
rate the water and to a small extent hy the secondaiy heat 
transfer due to 'the sensible heat contained in the water 
removed by evaporation. The heat transfer by evaporation 
is given by Hindi ey and Miner Lvl 

Qg = 2 h & (Tg) - e (la^J (2.21) 

where e (T^) is the partial vapour pressure at water sur- 
face temperature in mm of Hg, e partial vapour 

pressure at ambient temperature in mm of Hg and h is 
the heat transfer coefficient given hj equation (2.20). 
Also 

G = Gs ^^a^ ^ 

where e„ (T„) is the saturated vapour pressure at ambient 

s ^ 

tempera turo and R,H, is the relative himiidity at ambient 
temperature. 

Ttie rate of net heat loss, Q, from the surface 
is given by 


Q. 


0 . 


'■br 


"■ Q,. + 


Q 


( 2 . 22 ) 
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IJhen heat is being rejected by a power plant to a 
water body, the average surface temperature rises and the 
equilibrium condition is obtained when 

Qp = Qg + - Qg + or Q = Qp (2,23) 

where Qp is the rate of heat rejected by the power plant 
per unit area of the water body surface. 

2.3 One Dimensional Model 

Physically, a one dimensional model represents 
the circumstance where mixing in the water body in two 
directions is so vigorous that exchange of heat and mass 
need be considered only in the third direction. One 
dimensional horizontal model with energy exchange at the 
surface and variation in the longitudinal direction, is 
considered. Therefore, temperature in the vertical and 
lateral directions is assumed to be uniform. The theore- 
tical model, under steady state conditions, is developed 
below. 


2,3*1 Governing Equations 

The surface of the water body is assumed to 
consist of a number of finite divisions, each of width 
Ax. 'Ihe finite difference equations arc obtained by 
considering tlie energy balance for each element i.e., 

Energy input - Energjr output = Losses. 
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The governing equations for different elements 
have been developed in Appendix II. At intake and outfall, 
energy balance for half the elements is considered. The 
heat transfer due to conduction at intake and outfall has 
been assumed to be negligible. Thus three different 
equations have been developed, viz., one for the element 
at intake, one for the element at outfall and the third 
for all the intermediate elements. 

These eauations are solved alongwith the following 
boundary conditions, 

(a) Vr.lue of equilibrium temperature (i.e., temperature 
without heat rejection) throughout the water body is 
known, 

(b) Difference between the intake and outfall temperatures 
and the amount of heat to be rejected is knovm. This 
depends on the capacity of the thermal power plant and 
its efficj.ency. 

The solution is obtained by using the iterative 
method explained in section 2,1.2, 

2,4 Two Dimensional Model 

A two dimensional moded. represents the circum- 
stance where mixing in one direction is so vigorous 
that energy exchange need be considered only in the 



28 


other two directions, A two dimensional model with heat 
transfer and flow in the longitudinal and lateral direc- 
tions is considered# Thus a uniform temperature in the 
vertical diirection is assumed, with temperature variations 
in the longitudinal and lateral directions only. A theore- 
tical model, under steady state conditions, is developed 

DClOW, 


2.4.1 Governing Equations 

The surface of tire water bod;?" is divided into a 
number of square raeshes each of length Ax. The governing 
equations are obtained by considering the energy balance 
for each element as has been dj.scussed in the case of the 
one dimensional model. Appendix III outlines these 
equations. 

The equo/oions have been solved for the temperature 
field throughout the water body by using the iterative 
method. The velocity field obtained from the corresponding 
governing equations, given in Section 2,1, has been used 
in these equations. The following boundary conditions have 
been used. 

(a) Same as for one dimensional flo\r. 

(b) Same as for one dimensional flow. 

(c) The bouidaries of the water body are assumed to 

be adiabatic i.e., = 0 where n is a unit 

normal vector at the boundary walls. 
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The energ37' exchange between the water body and 
the environment was calculated after obtaining the tempe- 
rature field in the water body under steady state conditions. 
This energy was found to approximate to the heat rejection 
from the power plant. 
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RESULTS MD DISCUSSION 

This chapter presents the results obtained from 
the numerical solution of various governing equations 
developed in Chapter 2. Hie physical implications of 
these results have also been discussed. Of particular 
interest are the flow field, presented mainly in terms of 
stream function, and the temperature field, presented as 
isotherms. 

The two dimensional, horizontal flow patterns in 
the water body, under stead.y state conditions, has been 
obtained from the numerical solution of Equations (2,13), 
(2.1 5) 3-nd (2.l6), Figure 2 shows the streamlines in a 
water body, with inoake and outfall located on the same 
side, for the txTO limiting cases of inviscid and creeping 
flow as well as for flovfs at two other Reynolds number 
values, Ttie flow patterns for the inviscid and the creep- 
ing flow ci rcum stance v/i til intake and. outfall located 
on the ends of same side, are found to be sjmimetrical vxith 
respect to the centre line of the water body, whereas 
that for other Reynolds numbers is unsymmetrical due to 
the non3-inearity of the governing equations. A greater 
portion of the cooling water body is affected by the flowr 
in the case of inviscid flow as compared to other flows. 


Figure 3 shows the effect of aspect ratio of 
the water "body on the flow field for creeping flow. It 
is found that the portion of the water body, which is 
away from the outfall of the hot water, remains largely 
undisturbed. Hieref ore, in this case, tJrLs portion of 
the cooling water body does not take part significantly 
in the heat rejection process. It is found that the flow 
spreads out to only about half the distance between 
intake and outfall. Figure 4- shows the flow pattern for 
the creeping flow circumstance in a water body whose 
longitudinal dimension is much smaller than the distance 
between ulie intake and outfall locations. It is found 
that almost the entire water body gets disturbed in this 
case and, thus, helps in the rejection of more heat per 
unit area. The effect of locating the intake and outfall 
openings nonsymmetrically, on the flow pattern in a water 
body, is shown in Figure 5. 'fhe flow pattern becomes 
unsymmetrical even for the creeping flow. The portion 
of the water body, which does not fall between the intake 
and outfall locations, is found to remain largely undis- 
turbed and thus does not help in the removal of heat. 

Tho results obtained indicate circulatory, 
motion near the corners of tlie water body due to the 
formation of eddy currents. The strength of these cir- 
culatory motions is, however, found to be very weak and, 
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therefore, the same have nob been depicted on the 
Figures. 

Figure 6 shows the variation of velocity compo- 
nents u and V along the X - axis and Figure 7 shows 
their variations along the Y - axis, both for creeping 
flow, with intalre and outfall located on the same side 
of the water body. Figure 6 shows that the velocities 
have their maximum values near the outfall of the hot 
water and then decrease sharply away from it. The veloci- 
ties drop to almost ^ 0 % of the maximum value even before 
half the length is reached. Therefore, the portion of 
the water body near the outfall is very significant in 
the removal of the rejected heat. Figure 7 shows that 
the velocity component u decays to zero along the axis, 
midway between intake and outfall locations and then 
starts picking up the corresponding values symmetrically 
but in the opposite direction. 

The flow pattern in the water body with intake 
and outflow on the opposite sides of the water body is 
shown in Figure 8. Here also the inviscid flow affects 
a greater portion of the water body as compared to other 
cases. In this configuration, most of the water body is 
used, ilie effect of the aspect ratio of the water body 
on the stream function distribution for creeping flow is 
shown in Figure 9. It is found that the portion of the 


water body away from the outfall of the hot water is also 
disturbed and this results in more affective heat rejection. 
Figure 10 shows another configuration of intshe and outfall 
locations. It shows that the major portion of the water 
body remains largely undisturbed, if the intake and outfall 
are located on the centre of opposite sides. Hiis is, 
therefore, a very ineffective method of heat rejection to 
the water body. 

here, the results for the flow field, for flows 
with Reynolds number below 28.0, only, could be obtained. 
When higher values of the Reynolds number were used, the 
finite difference equations became unstable. For calcula- 
ting flow fields at higher Reynolds numbers, other methods 
such as the A.DI inobhod, employing a different finite 
difference scheiiio must be employed. In many of these 
methods, the transient problem is solved to reach the 
steady state solution. 

The numerical solution of the energy equations, 
to obtain the temperature field, is considered next. The 
net energy lost by the surface of the water body, Q, as a 
function of the surface temperature, wind velocity, cloudi- 
ness, relative humidity, etc., is computed by a numerical 
evaluation of the various components of energy exchange 
outlined in Section 2.2. The energy exchange has been 
analysed for a particular set of metaonolmjfginal parameters. 
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Figure 11 shows the effect of power plant capacity, 
i.e.j of bhe heat re;iccted per unit surface area of the 
water body, and bhe effect of eddy diffusivity on the intake 
temperature obtained from a one dimensional horizontal flow 
model. The intalie temperature increases with increase in 
Q as oxpeebed. Tlie average surface equilibrium temperature, 
t , which is defined as the surface temperature at which the 
net heat exchange from the surface is zero, is given by the 
inte 3 ?cept of the curve on the temperature axis. It is 
found thab the rate of increase in the intake temperature 
is small for lower values of Q. For higher values of Q, 
the intake temperature varies almost linearly. It is also 
found that the curve becomes slightly steeper with an 
increase in the value of eddy diffusivity. A Mgher value 
of eddy diffusivity results in greater energy transfer to 
the neighbouring fluid and would cause a larger rise in 
the intake texiiperature. Figure 12 shows the temperature 
decay in the direction of flow for various values of Q. 

As the hot water flows from the outfall to the intake, its 
tempera bure decreases due bo the energy lost to the environ- 
ment. For low values of Q, the temperature attains the 
equilibrium temperature before the flow reaches the intake, 
leaving the intake temperature unaltered. As the power 
level is increased, bhe increase in poorer being due to an 
increase in outfall temperature, the temperature throughout 
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the water body including the intake temperature increases. 
Figure 13 shows the effect of eddy diffusivity, or of 
turbulent mixing in the water body, on the temperature 
profiles. The intake temperature increases with an increase 
in the value of eddy diffusivity, as the increased value of 
eddy diffnsivity implies more mixing in the water body 
which increases the rate of heat transfer across the water 
body. Thus a large amount of thermal energy is transferred 
to thci adjacent layers and vhlch^ therefore, reflects a 
higher temperature at intake. 

The isotlienns in a two dimensional, horizontal, 
creeping flow are shown in Figure l4-. Uie temperature 
decays away from the outfall as expected, A tempeiature 
rise of about 3 at the intake is observed. The effect 
of outfall velocity on the tempera'ture field is shown in 
Figure The temperature difference, ^t, between 

outfall and Intake is taken as constant at 10 °C, As the 
power plant capacity, reflected in terms of higher velo- 
city, increases, the temperature at the outfall of the hot 
water also increases. Therefore, the isotherms shift away 
from the outfall in the water body for -ithl sc cincum stance. 

It is found, however, that the isotherms are more spread 
out as the power plant capacity increases. Thus the 
surface area of the water body significantly involved in 
the rejection of heat increases, Fig:'ure 16 shows the 
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effect of power plant capacity, i.e., of the heat rejected, 
on the inbake temperature. The rate of increase in the 
intake tempera bure, with respect to Q, is low for small 
values of (-j, as was found in the one dimensional flow. 

But for higher values oi' Q, the intal.o temperature varies 
almost linearly and at a hi.gher rate. The effect of 
Reynolds nuinher on the isothem distribution is shown in 
Figure 17. The isotheims have been shown i,d.th dimension- 
less temperature, T, which is defined as 



It is found that there is a slight shift of isotheims 
inwards indicating a higher rate of heat rejection with 
increased value of Reynolds number. Thus the intake tempe- 
rature decreases with an increase in Reynolds number, for 
the same amount of heat to bo rejected. Reynolds number 
may increase because of an increase in the outfall velocity, 
which must also have a lower At for given heat rejected and 
would give the trends observed, or a decrease in eddy 
diffusivity, which was seen to give a similar effect ear- 
lier. 

The results discussed in this chapter are predo- 
minantly for creeping flov;s as most of the flows we are 
concerned with in practice, for heat rejection to water 
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bodies, have low velocities vrlien averaged for the two 
dimensional model. The actaal problem is three dimen- 
sional. But to S'tud3r shallow water bodies or the effect 
of horizontal separation of intake and outfall, this 
model may be employed. 



CHAPTER h 


COHCLUSIOI'TS 

From tile discussions presented in Chapter 3, the 
results obtained are found to be physically reasonable. A 
flow configuration of considerable practj.cal significance 
is the one in wiiich the intake and outfall are located on 
the same side, with a horizontal separation. The flow 
pattern obtained for this configuration indicates the limi- 
tations posed on the effective use of the longitudinal 
dimension of the x;ater body, for heat rejection to the 
environment, mainly by the, distance between the intake and 
the outfall. The flow field for the configuration with 
intake and outfall on the opposite sides of the water body 
have also been studied and these show that the entire 
surface area of the water body can be used effectively, in 
this case, for heat rejection and, thus, resifLting in a 
lower ilsc in tl'c intake temperature, 

'ilie effect of heat rejection from a power plant 
is found to be a rise in the intake temperature, which 
shows a greater change for a high power plant capacity, 
as expected. The effect of an increase in the value of 
the eddy diffusivity results, as expected, in an increase 
in the value of intake temperature because the increased 
value of eddy diffusivity implies more turbulent mixing, 
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which leads to an increased rate of heat transfer across 
the water body. Ihe intake temperature is found to decrease 
marginall 3 r with an increase in Reynolds number over the 
range of low values of Re studied. This is expectedj 
since Re increases due to an increase in the velocity or 
in the outfall width or due to a decrease in the eddy 
dif fusivi by. 

The results obtained provide guidelines for limit- 
ing the power plant capacity for a given water body, with 
a particular horizontal distance available. The results 
are diroctl.y ap])licable to a shallow water body. Further, 
the results have been obtained for an idealized configura- 
tion, of vertical sides and flat horizontal base and top 
surface. To consider the arbitrary shape of the water 
body, the ' 1 % stroajiilinc may be obtained and a simplified 
rectangular geometry can then be taken in the proximity 
of this streamline. To account for an actual shape is 
neither econoMical nor advisable. In order to obtain the 
flow and temperature fields for a general three dimensional 
problem, which is usually the case in practice, it is 
reasonable to solve for the two dimensionaJ. vertical and 
the horizontal flow patterns, which will give results to 
adequate accuracy. The vertical flow pattern will indi- 
cate the depth that may be employed for the horizontal 
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model and the floxr spread of bhe latter gives the trans- 
verse ciijiiension for averaging in the former model. These 
results may also he employed in vertical one dimensional 
models for stucL3d.ng the transient and stratification 
effects. 


\ 
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APPENDIX I 


FINITE DIFFEEENGE EQUATIONS 

The finite difference eqtiations for two dimensional 
recirculation model are given helow. 

The equation (2.15) is represented in the finite 

difference form by 

( Tp +4' 4.iF .) 

¥ ___ i±U i-l. t ,1 ±1 l.U-l- 


hecomes 


The equation (2.16) in finite difference form 


8 A - B - C 
20.0 


where 
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® = ^\+1,3+1 +\+1,3-.1 +^-1,3+1 +^i-1,3-1 ^ 

and C = ( ?i^2,3 ^i-2,3 ^ ^i,0+2 ^ 

The finite difference form of equation (2.13) is 

given by 


where 
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APPENDIil II 


ONE DIMENSIONAL FLOW 


One dimensional flow, with temperature variation 
in the longitudinal direction only, is considered. 

The surfa.co area of bhe water body is divided into 
N number of finite elements, each, of length Ax and 
width L, The governing equations are obtained by consider- 
ing the energy balance for each element i.e., 

Energy Input - Energy Output = Losses 


The solution of energy balance equation gives 


^i ~ ^ "'e.O “ ^i+1^ “ 0^ 


where 

Qc 

DD 




AX 


U A 

if. v~ 


and tj_ represents the temperature at any point in the flow 
field except at intake and outfall. 


For finding the temperature at intake and out- 
fall, the energy balance is convSidered for half the element 
i.e., with length , and the conduction teim at intake 

and outfall is taken to be zero. Thus for each element, 
we get an algebraic equation as the governing equation. 
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APPENDIX III 


TWO DIMENSIONAL FLOW 


Two dimensional horizontal flow, with negligible 
temperature varic?tion in the vertica], plane, is considered. 

Hie surface area of the water body is divided into 
M parts, each ol length Ax, along the y axis and N 
parts, eacJi of length ax again, along the x - axis. The 
genera], energy equation for each element with two dimen- 
sional flow is given by 


3 t 3 1 

^ 3x 3y 


k_ ( t ^ t V 

V o O / 


pC. 


PCpH 


'p 3X^ 3y^ 

where Q denotes the energy transfer per unit area. 


The conduction term at intake and outfall are 
taken to be zero, in the above equation, while finding 
the intake and outfall temperatures. The above equation 
is written in the finite di.fference form for each element 
so as to give as many number of equations as are the number 
of elements. These equations are then solved in the same 
manner as was done for one dimensional flow. 
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inviscid Flow 


Creeping Flow 

Re =100 



FIG 2 FLOW PATTERN IN TWO DIMENSIONAL HORIZONTAL! 
FLOW FOR VARIOUS REYNOLDS NUMBERS , WITH! 

intake and outfall on the same side of 

THE WATER BODY 




y 


FIG. 3 FLOW PATTERN IN THE TWO DIMENSIONAL 

HORIZONTAL CREEPING FLOW, WITH INTAKE AND 
OUTFALL ON THE SAME SIDE FOR VARIOUS 
ASPECT RATIO'S OF THE WATER BODY 
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F(G.4 EFFECT OF REDUCING THE LONGITUDINAL 

DIMENSION OF THE WATER BODY ON THE 
FLOW PATTERN IN THE TWO DIMENSIONAL 
HORIZONTAL CREEPING FLOW; WITH INTAKE 
AND OUTFALL ON THE SAME SIDE 


LI ■ 






»Y 

FIG- 5 EFFECT OF LOCATING THE OUTFALL, AWAY FROM! 
THE EDGE OF THE WATER BODY, ON THE FLOW 
PATTERN IN THE TWO DIMENSIONAL HORIZONTAL 
CREEPING FLOW, WITH INTAKE AND OUTFALL ONi 
THE SAME SIDE 
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FIG- 6 VARIATION OF VELOCITY ALONG THE 

LONGITUDINAL DIRECTION FOR CREEPING 
FLOW, WITH INTAKE AND OUTFALL ON THE 
SAME SIDE 



/U. v/U 
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Velocity Component in X Direction — 
Velocity Component in Y Direction — 


FIG. 7 VARIATION OF VELOCITY ALONG THE 

TRANSVERSE DIRECTION FOR CREEPING FLOW 
WITH INTAKE AND OUTFALL ON THE SAME 
SIDE 



▼ - INVISCID FLOW 

V CREEPING FLOW 

FLOW AT Re=l0-0 


FiG 8 FLOW PATTERN IN TWO DIMENSIONAL 
HORIZONTAL FLOW , WITH 

ON THE OPPOSITE ENDS OF THE WATER BODY 



FIG. 9 FLOW PATTERN IN THE TWO DIMENSIONAL 
HORIZONTAL CREEPING FLOW , WITH INTAKE 
AND OUTFALL ON THE OPPOSITE ENDS, FOR 
VARIOUS ASPECT RATIOS 
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FIG. 10 THE TWO DIMENSIONAL HORIZONTAL CREEPING 
FLOW FOR INTAKE AND OUTFALL LOCATED IN 
THE MIDDLE OF OPPOSITE SIDES OF THE 
WATER BODY 



Intake Temperature 


58 


32-0 

31 0 


30-0 



FIG.II EFFECT OF POWER PLANT CAPACITY ON THE 

INTAKE TEMPERATURE FOR A ONE DIMENSIONAL 
HORIZONTAL FLOW MODEL 
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FIG. 14 ISOTHERMS IN THE TWO DIMENSIONAL 
HORIZONTA_L CREEPING FLOW AT 
Vh = 2-5 xlO^m^/sec 




FIG. 15 EFFECT OF A VARIATION IN OUTFALL VELOCITY 
ON THE ISOTHERMS IN THE TWO DIMENSIONAL 
HORIZONTAL CREEPING FLOW FOR AT=10°C 


Intake Temperature 
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FIG. 16 EFFECT OF POWER PLANT CAPACITY ON 
THE INTAKE TEMPERATURE IN THE 
TWO DIMENSIONAL HORIZONTAL CREEPING 
FLOW AT Vh = 2-5 xIO'^ m^/sec 



FIG. 17 EFFECT OF REYNOLDS NUMBER ON THE 

ISOTHERMS IN THE TWO DIMENSIONAL HORIZONTAL 
FLOW AT VH = 10^m^/sec 




